Quantum reflection: The invisible quantum barrier 
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We construct the invisible quantum barrier which represents the phenomenon of 
quantum reflection using the available data. We use the Abel equation to invert the 
data. The resulting invisible quantum barrier is double- valued in both axes. We 
study this invisible barrier in the case of atom and Bose-Einstein Condensate reflec- 
tion from a solid silicon surface. A time-dependent, one-spatial dimension Gross- 
Pitaevskii equation is solved for the BEC case. We found that the BEC behaves 
very similarly to the single atom except for size effects, which manifest themselves 
in a maximum in the reflectivity at small distances from the wall. The effect of the 
atom-atom interaction on the BEC reflection and correspondingly on the invisible 
barrier is found to be appreciable at low velocities and comparable to the finite size 
effect. The trapping of ultracold atom or BEC between two walls is discussed. 

^ 1 I. INTRODUCTION 

Quantum reflection is an intriguing phenomenon which has been under experimental 
scrutiny in the last few years. Several papers have been written on the subject. Back 
in 1983, V. U. Nayak, D. O. Edwards and N. Masuhara [l| measured the scattering of 
■^j- ■ 4 He atoms grazing the liquid- 4 He surface and found that the reflectivity approaches unity 
as the velocity tends to zero. Berkhout et al. 0] measured the quantum reflection of 
H atoms from a concave spherical mirror. Yu et al. [1] revisited the experiment of [l[ 
and measured the sticking probability. More recently, [4| looked at the phenomenon of 
quantum reflection of very slow metastable neon atoms from a solid silicon surface and 
a BK7 glass surface. More data were collected in the following years. Sodium Bose- 
Einstein condensate (BEC) containing 3 x 10 5 atoms reflected at normal incidence from 
a silicon surface at incident velocities of 1 - 8 mm/s were reported by 0. Extension 
of this study to quantum reflection of a 23 Na BEC containing 10 6 atoms at velocities 
below 2.5 mm/s on a dilute silicon surface demonstrated the effect of the atom-atom 
interaction and the extended density profile as the reflectivity saturated at 0.6 rather 
than 1.0 [fjj. The observed reflectivity is explained as quantum reflection caused by the 
attractive Casimir-Polder (CP) atom-wall potentialM] The significant reduction in the 
reflectivity of BEC at very low velocities allegedly [8j arises from the modulation in the 
BEC density profile which ensues owing to the formation of standing waves resulting 
from the superposition of incident and reflected matter waves. As a result, dynamical 
excitations, such as solitons and vortex rings, are created in the BEC which fragment it 
and disperses its atoms. Such work on quantum reflection of BEC's is important for the 
understanding of the Bose-Einstein condensate stability near surfaces [qJ. Several recent 
theoretical works on the general foundation of quantum reflection have been published 
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The observed reflectivity is explained as quantum reflection caused by the attractive 
Casimir- Polder potential [7J. The reflection is accordingly classically forbidden. The 
Wentzel-Kramers-Brillouin (WKB) approximation would not work in such a case since 
the variation of the local de Broglie wavelength is comparable to itself, regardless 
to the sign of the interaction. This, however, does not preclude the use of adequate 
perturbation methods to calculate very approximately the reflectivity [Io|. On the other 
hand, reflection implies tunneling and one can thus invert the data on reflectivity, R(E) 
(or tunneling, T(E) = 1 - R(E)) using WKB type formulae [la ], to obtain the width of 
what might be called the invisible quantum barrier. Such a procedure would avoid the 
use of ill-defined concepts such as the quantum potential jioj ] . The purpose of this paper 
is to use the inversion procedure to obtain the invisible barrier responsible for quantum 
reflection. We apply this to the available data on quantum reflection of atoms [4J and of 
BEC [5|, |6[ . We also consider the resulting invisible barrier in the context of the trapping 
of BEC between two walls. 

The paper is organized as follows. In section II, the phenomenon of Quantum Re- 
flection is reviewed and the Casimir-Polder potential is discussed. In Section III, the 
Invisible Barrier is derived using the Abel equation for existing data on atom-wall and 
BEC-wall reflection. In Section IV, the Invisible Barrier is obtained for a BEC using 
a one dimensional Gross-Pitaevskii equation. The free space trapping of BEC between 
two walls is also considered. The trapping occurs within the confines of the two invisible 
barriers which sit at an appreciable distance far from the walls that create them. Finally, 
in Section V, several concluding remarks are made. 



II. QUANTUM REFLECTION 

The phenomenon of quantum reflection occurs whenever the local de Broglie wave 
length of the moving particle is comparable to the distance over which the potential 
varies rapidly, regardless to whether the potential is repulsive or attractive. For an 
attractive potential, such as the one between an atom and a wall ( the Casimir-Polder 
interaction), the quantum reflection is a classically forbidden process, just like tunneling 
through a barrier. In this section we give a short account of Quantum Reflection, in 
connection with the CP interaction. 

The CP, atom- wall, interaction, which has the same physical vacuum fluctuation origin 
as the famous Casimir force between parallel walls separated by a distance L, , Vc{r) = 
1 is derived in 0] and can be written as an integral ( see, e.g. 14 . [TBI] ) . 
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where a(iu) is the dynamic electric dipole polarizability of the atom evaluated at the 
imaginary frequency iu and a/ s is the fine-structure constant. The quantity a{iuj) can 
be evaluated using sum rules and it can be written as, 

a{iu) = S n — fi" 2 (2) 

(E n - E Q y + u) 1 



where the sum S n runs over the discrete dipole states and the continuum scattering 
states. In the above, the dipole oscillator strength f n is given by 
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In the above equation N represents the number of electrons. Further, the quantity 
S n fn is just the dipole sum rule wiV, fl~5| . 

The asymptotic form of the CP potential is invariably written as, 
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where C4 has the value,e.e;., 9.1 x 10 -56 J.m 4 for sodium atoms incident on a silicon 
surface and A = 590 nm, [5|, [6( 

The result of the Quantum Reflection measurements cited above show a strikingly 
robust quantum reflection at very low velocities. The presence of such quantum reflec- 
tion raises the question about the transmission to the solid wall, where the atoms are 
adsorbed. The transmission, or tunneling implies the presence of an invisible quantum 
barrier "sitting" at a distance close to the wall. Can one determine this barrier from the 
available data on quantum reflection? How does a BEC behave compared to a single 
atom? The answer to these queries is the thrust of the present paper. 

Knowing the reflection coefficient, the transmission one is simply the compliment to 
unity, viz T(E) = 1 - R{E). 

In the following we use the generalized WKB form for the tunneling probability, 
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Above, ri and r 2 are the inner and outer turning points which are the roots of the 
equation V(r) = E. Further, the above form of the tunneling probability is exactly 1/2 
at the top of the barrier as the exact solution of the Schrodinger equation so requires. 

Knowing T(E) from the data, one can invert it to obtain the barrier. This procedure 
is intimately related to the solution of Abel's theorem and the corresponding classical 
mechanics problem [171 ] . The equation that does this is given by 
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where B is the height of the barrier, generally defined through T(B) = R(B) = ~, which 
results from JJj 2 dr^JV{r) — B = 0. This last relation defines theoretically the height of 
the barrier, while the former supplies an "experimental" definition of this height. 

The above equation, which we shall refer to as Abel formula, has been used in a 
variety of forms to investigate inversion in the WKB type description of bound systems 
by [l8(. It has also been used to extract the barrier responsible for tunneling into the 
classically forbidden region in the sub-barrier fusion of nuclei [2(| . The resulting potential 
which results from the inversion of the data is multivalued function. Of course, by its 



nature, Eq. (6) dictates that the thickness function, r 2 — r 1; is determined by R(E), or 
equivalently, T(E) for V < E < B, namely in the energy region where R(E) > 1/2. 



III. INVISIBLE QUANTUM BARRIER 



Cole and Good 18]have suggested using the Abel formula to invert the semiclassical 
reflectivity and tunneling probability to obtain the thickness of the barrier in cases of 
reflection in the classically allowed region. Here we have used the Abel formula to obtain 
what we may call the invisible quantum barrier (INV) responsible for the quantum 
reflection (QR), reflection in the classically forbidden region, of single sodium atoms [H 
from a silicon surface, and for the QR of a Bose-Einstein condensate composed of a 3 x 
10 5 sodium atoms in a trap [E| and of a million such atoms 0]. Although the interaction 
between BEC and a surface is less well understood than that of a single atom with a 
surface, we take the practical view of using the same type of the CP interaction as was 
done in jl, @, When using Eq. (6), we have fitted the data points with an analytical 
function of the velocity and used this function in the inversion. Of course the inversion 
alluded to above would only supply the barrier value as a function of its thickness, 
Vi nv (xi — X2). Of course this procedure is applicable at energies where the reflectivity is 
larger than |. 

For energies greater than the height of the barrier, R(E) < |, one can use the 
parabolic barrier approximation valid near the position, Rb, of the height, B, of the 
barrier, V inv (r) « B - \\^^\ r =R B {r ~ r b) 2 , which gives R(E) = 1+cxp[ l (E _ B)] , where 

hu measures the curvature of the assumed inverted parabolic barrier, 
where m is the mass of the atom. This procedure is quite common in the theory of 
nuclear fission By fitting the data points for R(E) < | we easily obtain the values 
of the height ( which can also be read off from the condition R(E) = | ) and, more 
importantly, the curvature of the invisible barrier. The height and curvature parameter 
of the barrier for the data on Ne atoms reflection off a silican wall, [4j are, respectively, 
B = IMnK and ^^| r=i?s = 0.011 w ?u m i , and on a BK7 glass wall, 8, B = 0.05nK 
and d JJf^ \ r =R B = 0.015 10 -^ m2 ■ In the case of the Na BEC reflection data of [6] we find 

for these parameters the values, B = 23.7nK and d JJf^ \r=R B = 17.0 10 -u m2 ■ Clearly the 
barrier for BEC reflection is more than an order of magnitude higher and thinner than 
for single atom reflection. 

The Abel formula was employed using the available data for R(E) > \ and 
extrapolating these to lower energies using an appropriate fitting function. The fit 
formula which we have employed is of the general form R(E) = exp [— CE D ]. Such a 
functional form reproduces the over all trend of the data of [4] and to a large extent 
those of the BEC reflectivity of [6]. In the case of reflection of neon atoms from a 
silicon wall of [4|, we find C = 0.81 and D = 0.71, whereas for reflection of neon 
atoms from BK7 glass wall, C = 0.92, D = 0.29. In the case of the sodium BEC 
reflection from a silicon wall 0, we obtained C = 0.23 and D = 0.39. This latter case 
could account for the small reflectivity (higher velocities) and misses altogether the 
observerd saturation at smaller velocities. The inversion alluded to above would only 



supply the barrier value as a function of its thickness, V inv (xi — x 2 ). Such potentials 
are shown in Fig. [TJ The barrier heights and curvatures obtained from the inversion 
procedure are: B = 1.03nK and d JJf'* \r=R B = 0-0086 1Q -S m 2 for neon on silicon wall,j3], 
B = 0.0A8nK and ^P-\ r =R B = 0.0108 10 ?£ m2 for neon on BK7 glass wall, fl, and 

B = 23.71nK, d H r - ) | r= R H = 17.64 irv lffi 2 fo r sodium BEC reflection from silicon wall 
[6j. The values of the curvatures are quite consistent with those obtained using the 
parabolic approximation above. 

The result of the calculation of the BEC Vi nv using the time-dependent Gross- 
Pitaevskii equation is given below. Before we discuss these results, it would be useful to 
obtain the invisible barrier as a function of the distance between the reflected entity and 
the silicon wall. For this purpose we need another equation involving the two turning 
points, Xi and x 2 . One possible suggestion is to relate their product to the square of 
the virtual turning point, r(E), related to the Casimir- Polder potential Vcp, namely 
E = Vcp{r(E)), and r(E) 2 = x\ • x%. Through this relation we get the real turning 
points x\ and £ 2 separately, and accordingly we can construct the invisible quantum 
barrier as a function of the distance from the wall. This procedure is not unique as a 
different relation involving the turning points and the Casimir-Polder potential would 
yield a different shape for the invisible barrier. Not having available another plausible 
relation, we proceed and use the product relation above. The result of such a calculation 
is presented in Fig [2J We see clearly that the barrier is a double valued function in 
both axes. This feature may be a reflection of the use of an energy- dependent relation 
between the turning points and r(E). It would be certainly of value to explore other 
relations, such as X\ + x<i = 2 • r(E). 



IV. INVISIBLE QUANTUM BARRIER FOR BEC 



We now turn to the study of quantum reflection of an one dimensional Bose-Einstein 
condensate from the silicon surface. Several theoretical studies of ID BEC have been 

]). It has been argued that by lowering the dimension to ID, 



published ( see, e.g. [22 



the transition temperature for BEC formation is increased. These systems are subtle and 
exhibit features not encountered in 2D or 3D. Bearing in mind the intrinsic differences 
between the ID and 3D BEC's, we shall, nevertheless, discuss the quantum reflection of 
ID BEC for the purpose of simplicity. Would the condensate suffer reflection just as a 
single atom? In a way, to answer this question is in line with a broader one connected with 
the quantum behaviour of the motion of mesosco pic systems in general. One is reminded 
here with the pioneering work of M. Arndt 



24. 25 



on interference pattern of Cqq and 
other heavy molecules as they pass through a grating. As in [g], we consider N 23 Na atoms 
Bose-Einstein condensed into the ground state in an anisotropic trap, Q p ^> Q x , where 
Q p is the trap frequency in the transverse directions and Q x the frequency in x-direction. 
Both the quantum reflection and the interaction between atoms will not excite the motion 
in the transverse direction. With these assumptions, the condensate is described by the 
one-dimensional Gross-Pitaevskii (GP) equation 1261 1 . In fact one dimensional BEC has 
been formed and studied in 0, [IT 
2^ 



The ID GP equation we shall solve here is 
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where g = N x g 1D with N the number of atoms and gin = 2hQ p a, where a is the 



s-wave scattering length, which can be tuned through the Feshbach resonance [27J. The 



condition of validity of the above formula for gin is: a <C \/~jr 2fl- Notice that the 



number of atoms appears in the strength of the interaction g as the wave function has 
been normalized to unity, 

\tp\ 2 dx = 1. (8) 

The external potentials include the harmonic trap in x-direction and Casimir-Polder 
potential induced by a silicon surface (1, 0], V(x) = V tr ap + V cp . The harmonic trap is 
Vtrap = fiVL 2 x x 2 /2. In this paper, the trap frequency is set to be Q x = 2n x 3.5Hz, in 
the range of values of recent experiments 0, 0, 2$. The Casimir-Polder potential has 



the form of Eq. (4) . Similar to [8J, the model Casimir-Polder potential is separated 
into two regions. If x < d — 5, we have the normal Casimir-Polder potential. Here 
5 = 0.15/zm is a small offset to silicon surface. When x > d—5, we use a complex potential, 
K — V cp (5) — i(x — S)Vi modeling adsorption at the surface with Vi = 1.6 x 10~ 26 J- m^ 1 . 

Recently experiments on the reflection of condensates were realized in three dimen- 
sional traps with number of condensed atoms reaching 10 6 [H, 0]. However, the number 
of atoms which can be condensed in one dimensional trap is reduced due to thermal 
and quantum fluctuations. We use the number of atoms in the BEC of a recent exper- 
iment |29|], Nia = 1.5 x 10 4 . At first we "prepare" a ground state of condensate in the 
harmonic trap. The center of trap was shifted to the Si surface suddenly by a distance 
As. Therefore the condensate is suddenly put at a high potential: AV = MQ 2 Ax 2 /2. 
This induces the condensate to travel towards the surface, see Fig. [3J At the surface, the 
incident velocity is approximated by v x ~ fl x Ax. We use a fourth-order Runge-Kutta 
method in the interaction picture to perform the time propagation of the solution of the 
Gross-Pitaevskii equation. After one complete reflection, we define reflection probability 
as the number of left-traveling atoms divided by the number of incident atoms. Different 
from the single atom reflected by the surface, the condensate has a finite spatial exten- 
sion. This yields a minimal but nonzero velocity of incident condensate. In addition, 
during the reflection process, atoms are lost continuously. The nonlinear term in the 
Gross-Pitaevskii equation is also reduced accordingly since it is proportional to the total 
number of atoms in the condensate. 

The simulation is presented in Fig. HI When the velocity is reduced, the reflection 
probability approaches unity. However, due to the facts mentioned above, the velocity 
can not reach zero. So, in principle, we get a less-than unity maximum reflection, similar 
to the experimental results of [6|. With the reflectivity calculated, we can now obtain 
the invisible barrier following the procedure detailed above. 

We use the Eq.(6) to get the barrier at different thicknesses. The corresponding 
Vi n Jj) is shown in Fig. [HI and it resembles very well the one obtained from the data 
of [6j, Fig. [2J This attests to the consistency of our calculation. More importantly, 
it shows clearly that a BEC behaves very much like a single atom, albeit with a finite 
size. The size effect is felt when the BEC comes very close to the wall, as has already 
been mentioned. We have also investigated the effect of the strength of the atom-atom 



interaction, g, on the reflectivity of the BEC. In the numerical simulation, we scale the 
energy and length with ground state harmonic oscillator energy, E = hQ x and harmonic 
length Z = \jTi/mVl x . In this way, we solve the equation with scaled interaction strength 
g s = 0.0 ( ideal gas ), g s = 50, and g s = 200.00. The resulting reflectivity and the 
corresponding Vi nv {r) shows appreciable sensitivity to the value of g, and accordingly on 
the scattering length, a, Fig. [6l which is comparable to the finite size effect alluded to 
above. We have also calculated the parameters of the barrier using the inverted parabola 
approximation of Hill- Wheeler. We find, for the height B = lllnK and the curvature 
d \ r =R B = 195(g = 0.0), 141(g = 200) 10 -u m z ■ These values are to be compared to 
the ones extracted from the data, mentioned above. The calculated ones are a factor 5-6 
larger than the ones obtained from the data. This discrepancy could be traced to the 
several limitations inherent in our ID GP simulation. 

It has been recently suggested that that one could trap cold atoms and BEC by 



quantum reflection from two walls [33, |34| . Our findings above give a specific mechanism 
for such trapping. The cold atom or BEC will be in the confines of the two invisible 
barriers that such walls will generate. As long as the velocity of the atom or the BEC 
does not exceed the top of the invisible barrier, they will be trapped in free space without 
touching the walls in accordance with the conclusions of 33|, 3J] . The wall-atom- wall CP 
potential can be written as [15| 
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where z is the distance of the atom from the mid-point between the walls, taken to be 
separated by L. The profile of the above WAW potential is shown in Fig. [5j For the 
purpose of completeness we calculate the corresponding double invisible barrier when the 
BEC is placed between two walls. 

A schematic picture of the confining double-barrier in the case of the sodium BEC is 
shown in Fig. [61 One sees clearly that the double-barrier defines the confining space to be 
the interval ( -22, +22) x 10~ 7 m, while the actual walls are situated at (-26, +26) x 10~ 7 
m. The space between the barrier and the wall is completely inaccessible. This, however, 
depends on the energy, as the invisible barrier is energy-dependent. In the case of sodium 
BEC |6[, the extracted critical velocity which the BEC must have in order to reach this 
forbidden region corresponds to an energy of about 24.0nK. Our ID GP calculation 
gave a height which is 5 times greater. This is due to the inherent limitation of the ID 
simulations. We should also warn the reader that our confining box-like potential was 
constructed with using the plausible condition on the turning points, x\ x xi = r(E) 2 , 
with r(e) being the virtual turning point associated with the Casimir-Polder potential. 
Another condition would give a different confining region. 



V. CONCLUSIONS 

In conclusion we have inverted the recent data on reflectivity of neon atoms and of 
sodium Bose-Einstein condensates to obtain the invisible quantum barrier. We have 
found that this new object is double valued in both axes. The height and curvature of 
the invisible barrier were determined from the data. The value of the height we obtain 
form the data of [6] is B = 2AnK, while for the neon atoms of [4j it is l.OnK. The former 



is high enough that makes the trapping of BEC quite feasible at energies smaller than B. 
Such trapping could be potentially of great value in the fabrication of atom laser devises, 
as the BEC will maintain its integrity without suffering from absorption-adsorption at 
the walls. Further, we have verified that a BEC, though mesoscopic and contains millions 
of atoms, behaves, to a large extent, like a single, albeit large, atom when reflected from 
a solid wall. This finding corroborates those of 0] on the persistence of quantum 
behavior of large mesoscopic objects and is in line with the results of the recent grating 



diffraction experiments on buckyballs and other large molecules [2J, |25 
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FIG. 1: The invisible barrier vs. its thickness for the data of [4|, |a|. The experimental re- 
flectivities were adjusted with a simple function and used in Eq. (6). Insets: Fig la., The 
reflectivity data versus the normal incident velocity on the Si(l,0,0) surface m. Fig lb., The 
reflectivity data versus the normal incident velocity on the BK7 glass surface, [4]. Fig la, Re- 
flectivity of sodium BEC data versus incident velocity. Data points, 0], correspond to a sodium 
BEC confined in a magnetic trap with trap frequencies 27r(2.0, 2.5, 8.2) Hz (full squares) and 
27r(4.2, 5.0, 8.2) Hz (open triangles). See text for details. 
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FIG. 2: The invisible barrier vs. the distance, x, for Ne reflection from a silicon wall ( at 
x = 0.0 ), from a BK7 glass wall, [4] and for Na BEC reflection from a silicon wall 0,0. The 
dashed curve indicates the Casimir-Polder potential. See text for details. 




FIG. 3: (Color online.) A schematic diagram of the potentials. The blue curve is the harmonic 
trap. A ground state condensate was created in the trap. The red curve is the Casimir-Polder 
poential. In the numerical calculation, we set a finite offset to the surface(see text). The black 
curve indicates the imagnary part of the absorbing/adsorbing potential. If atoms reach this 
region, they will be adsorbed at the surface or inelastically scattered away from the trap. 
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FIG. 4: (Color online.) Quantum reflection of one dimensional sodium Bose-Einstein condensate 
from a silicon surface. Velocity is measured in mm/ s. Three values of the atom-atom interaction 
strength were considered: g s = 0.0 (dashed line), g s = 50.0 (dashed dotted line) and g s = 200.0 
(solid line). See text for details. 
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FIG. 5: The profile of the wall-atom-wall potential according to Eq. (14). The walls are situated 
at ( -26, +26)x 10~ 7 m 




FIG. 6: (Color online.) Top: The invisible barrier of a silicon surface ( at x = 0.0 ) for a sodium 
Bose condensate for two values of the interaction strength, g s = 0.0 (dashed) and 200.00 (full). 
Here we use the units: distance in meter and energy in nK. Bottom: The same as the above 
figure but with two silicon walls situated at x = (—26, +26) x 10~ 7 m, schematically showing 



